Molecular Ecology

Trait-selective harvesting by fisheries can impose strong selective pressures on fish populations, driving changes in life history traits affecting fisheries productivity and ecosystem functioning. While the genetic consequences of harvesting have been extensively studied, the extent to which phenotypic variation reflects genomic evolution versus environmentally-induced plasticity remains unclear. Epigenetic mechanisms, such as DNA methylation, may mediate between these processes, serving as a rapid and reversible response to the selective pressures imposed by harvesting. In this study, we implemented an improved laboratory and bioinformatics protocol, epiGBS3, to examine genomic variation and DNA methylation patterns in the marine fish Xyrichtys novacula. The study spanned three replicated geographical areas each comprising two adjacent locations: an intensively exploited fishery and a no-take Marine Protected Area (ntMPA). A nested analysis design across the three areas revealed strong gene flow and no evidence of genetic structure. Nevertheless, nucleotide diversity was significantly reduced in fisheries relative to ntMPAs. We also found that DNA methylation levels differed between protected and exploited sites after controlling for age, suggesting that fishing may influence epigenetic changes independently of fisheries-induced age-truncation effects. This represents one of the first lines of evidence that fisheries can potentially shape epigenetic variation, supporting DNA methylation as contributor to local adaptation under high gene flow and strong anthropogenic selection.

Introgressive hybridization between wild and domestic animals is a widespread phenomenon with important implications for genetic diversity, local adaptation, and conservation management. The causes and consequences of this process are poorly understood. In Australia, hybridization between dingoes and domestic dogs presents a dual conservation challenge, threatening the genetic integrity of dingoes while allowing potential adaptive introgression. To investigate the environmental drivers of this process, we analyzed high-density SNP array data in 390 dingoes and 396 domestic dogs. Dingo populations showed regional genetic structure and were clearly differentiated from domestic dogs. Using local ancestry inference and genome-environment association analyses, we found low levels of dog introgression in dingoes from remote areas in Central and Western Australia, and moderate levels in Eastern and Southern populations. Climatic variables (maximum temperature of the warmest month, mean temperature of the driest quarter) and the Human Footprint Index (reflecting density of human populations and environmental modifications) were significant predictors of introgression. We identified four genomic regions with overrepresented dog ancestry, including a large introgressed block on chromosome 27, which contained an olfactory receptor gene showing signatures of positive selection, suggesting adaptive introgression. In addition, a chromosomal inversion previously described in dogs and absent in dingoes was initially identified as an introgressed block. We also detected eight genomic regions nearly free of dog ancestry, suggesting purifying selection against maladaptive variants. Together, these results highlight the complex interplay between introgression, human influence, and local adaptation in dingoes, offering valuable insights for conserving the evolutionary potential of this apex predator in increasingly modified landscapes.

Kelp forests are widely distributed along temperate and polar coastlines worldwide and are among the worlds most productive and diverse marine ecosystems. Yet, due in part to ocean warming, they are declining and even disappearing in many parts of the world. While genomic tools can identify local adaptation and predict species responses to global change, these predictions have rarely been validated in the field, hampering their widespread use in conservation practice. Here, we applied a seascape genomics approach to investigate environmental adaptation in the two main canopy-forming species of the Northeast Pacific, Macrocystis tenuifolia and Nereocystis luetkeana. We leveraged whole-genome sequences of 598 individuals across 94 sites along the British Columbia and Washington coasts, together with 37 environmental variables. Both species showed genomic signatures of local adaptation, with distinct environmental drivers shaping adaptation in each species despite their co-occurrence across much of the studied area. Using gradient forests, we modelled the genetic turnover across environmental gradients and predicted populations vulnerability (genomic offset) under projected environmental conditions. Genomic offsets differed greatly among regions and were positively correlated with kelp declines observed to date, especially in Macrocystis, validating the link between genomic models and outcomes in the field and allowing us to translate genomic predictions into an ecologically meaningful metric: the risk of extirpation under global change. Our models predict that assisted migration could significantly attenuate kelps vulnerability to global change. Across environmentally heterogenous coastlines, short-distance migration can often substantially reduce future genomic-environmental mismatches, but in many cases, long-distance migration would be most beneficial. Our results highlight the potential of seascape genomics to predict vulnerability of populations to global change. Importantly, the validated link between our genomic models and ecological outcomes allows quantification of climate-driven extirpation risk and can inform conservation strategies to improve the resilience and sustainable management of these vulnerable ecosystems.

Codiversification often arises when hosts and their endosymbionts share a linked evolutionary history, exhibit vertical transmission, or share ecological and biogeographic processes. Most studies on the codiversification of carpenter ants (genus Camponotus) have focused on the co-phylogeny of hosts and endosymbionts across multiple species; however, no studies have examined the intraspecific population-level phylogeographic patterns of codiversification within Camponotus. California has been a geographic focus for phylogeographic studies due to its high endemism and complex geographic structure, and Camponotus laevigatus is a carpenter ant primarily found there. Here, we used whole-genome sequencing from C. laevigatus and its endosymbiont, Blochmaniella to investigate phylogeographic patterns of host-endosymbiont codiversification and estimated kinship of ants sampled near one another. We identified three phylogeographic clusters and isolation-by-distance analyses indicated a positive relationship between genetic and geographic distance in C. laevigatus and Blochmaniella. Using estimates of effective migration surfaces, we found that the Central Valley in California acts as a significant barrier to gene flow among populations. Our phylogenetic analyses revealed the congruent phylogenies of C. laevigatus and Blochmaniella, supporting codiversification. We also estimated kinship among individuals from the same and nearby sampling sites; kinship results indicated full-sister relationships among individuals from the same sampling site, except for three pairwise comparisons, and foragers from nearby sampling sites displayed some shared kinship. Lastly, our demographic analysis revealed a Pleistocene divergence, highlighting the role of Quaternary climatic cycles in shaping the population structure of C. laevigatus.

Population genetic theory predicts that natural selection will be more efficient in large than small populations because genetic drift reduces the efficiency of selection in small populations. Small populations adapting to new environments can also be expected to evolve higher recombination rates to facilitate adaptation as well as to dissociate and purge harmful mutations. We tested these hypotheses (1) by investigating differences in the strength of association between nucleotide diversity ({pi}) and recombination rate across the genomes of nine-spined sticklebacks (Pungitius pungitius) from four small freshwater (mean Ne {approx} 2 578) and four large marine (mean Ne = 86 742) populations, as well as (2) by comparing recombination rates between small and large populations using population specific linkage maps. We found the predicted positive correlation of{pi} with recombination rate from all but the smallest freshwater populations, suggesting prevalent linked selection even after accounting for variation in GC/CpG content, and gene density. Mean recombination rates did not differ between freshwater and marine populations, except that the smallest Ne freshwater population exhibited significantly elevated recombination rate. GWAS analyses suggested a polygenic basis for recombination rates. These results suggest an important role for linked selection in reducing{pi} in low recombination regions especially in large populations. Moreover, as predicted by theory, at least one of the small freshwater populations appears to have evolved a higher recombination rate than its marine ancestors.

Restoration and management of natural populations often assume that local genotypes are best suited for transplantation to their local environment. Prioritizing a single local genotype, however, contrasts with the framework of maximizing intraspecific diversity to increase population resilience to environmental change. Local populations may also become maladapted to a rapidly changing environment, motivating alternative frameworks that instead minimize environmental distance between source and transplantation sites. Here, we tested the predictive power of the local is best, maximize intraspecific diversity, and minimize environmental distance frameworks on the survival and growth of Eastern oyster (Crassostrea virginica) genotypes in field common gardens that differed in salinity and disease pressure. Although a genome scan revealed patterns of adaptation to disease, heat stress, and salinity among source populations, we did not find strong support for the local is best framework: geographically distant southern genotypes performed comparably to local selection lines and a local wild population. Higher genetic diversity within monocultures was associated with higher survival, yet highly diverse polycultures survived at lower rates than the best-performing monocultures, providing mixed support for the maximize intraspecific diversity framework. Temperature and salinity of the environments-of-origin of parents predicted the survival of their offspring in common gardens, but the relationship between survival and environmental distance was context-dependent, leading to mixed support for the minimize environmental distance framework. Together, these results demonstrate that no single framework reliably predicted transplantation success, suggesting that effective management strategies may need to integrate genomic and environmental lines of evidence to guide genotype selection.

Heat Shock Protein 90 (HSP90) functions as an evolutionary capacitor, allowing populations to store cryptic genetic variation that can be released under stress. While former studies have described the release of morphological variation, its behavioural consequences remain unexplored. In the red flour beetle, Tribolium castaneum, HSP90 inhibition released a phenotype with much smaller, less defined eyes that confers fitness benefits in continuous light and was subsequently assimilated. We hypothesized that altered eye morphology affects light perception and thereby changes light-dependent behaviours. To test whether phenotypes released via evolutionary capacitance can beneficially alter behaviour, we examined locomotor activity rhythm entrainment to light-dark cycles as well as individual and group light choice behaviour. Males of the reduced-eye phenotype exhibited a diminished startle response to sudden light exposure in locomotor activity assays. We also found reduced negative phototaxis in groups of beetles with reduced eyes. This modified behaviour, indicating reduced light sensitivity, may stem from impaired light perception caused by altered eye morphology. Lower light sensitivity could be beneficial under stressful environmental conditions by promoting the exploration of alternative niches. Therefore, this study provides the first evidence for potentially beneficial behavioural changes in a HSP90-released phenotype, reinforcing HSP90s role as an evolutionary capacitor.

AO_SCPLOWBSTRACTC_SCPLOWRockfishes (genus Sebastes) show extreme variation in longevity among closely related species, but the evolutionary history of this young radiation is highly complex. To unpack these relationships and to associate genotypes with phenotypes, we quantified genealogical discordance among 55 Sebastes species and implemented a phyloGWAS framework that incorporates discordant gene histories into genotype-longevity association tests. We found that genealogical discordance is extremely high: the inferred species tree topology differed among several ILS-aware methods, with most internal branches having low concordance factors regardless of which method was used. Nevertheless, some phylogenetic structure was shared by all inferred species trees. We used simulations to assess the statistical properties of phyloGWAS applied to complex traits using different genetic relatedness matrices (GRMs) and under varying levels of discordance. Adding an accurate GRM reduced false positives relative to a model without relatedness, but GRMs only modestly increased power to detect true positives. Using multiple approaches on the Sebastes data, phyloGWAS identified several variants associated with longevity. Our results indicate that extreme genealogical discordance is a core feature of Sebastes evolution and that phyloGWAS can help in connecting genotype to phenotype under these conditions.

Identifying spatial and temporal patterns of connectivity among populations is fundamental to marine ecology, evolutionary biology, and fisheries management. Yet, due to large population sizes and low genetic differentiation among populations, empirical quantification of population connectivity across a species entire range has not been achieved for an open-coast marine organism. Here, we leverage experimental transcriptomics to develop a genotyping-in-thousands by sequencing (GT-seq) panel to support assignment of recruits of the kelp forest gastropod, Kellets whelk (Kelletia kelletii), collected across the species biogeographic range. Over a three-year period, we identified high self-recruitment in the historical range (100%) and low self-recruitment in the expanded range (10.53 - 13.73%). Additionally, self-recruitment within the expanded range generally increased with recruit age, from 27.14% at 0.93 years to 43.40% at 1.93 years, indicating that the locally spawned individuals were more likely to survive to older ages than migrants from the historical range. Together, these results reveal limited self-recruitment in the expanded range and suggest that a post-settlement selective filter contributes to differential survival in a high gene flow marine system.

Successful establishment of a species in a new range is a useful way to understand the impact of demography and selection on the evolution of globally distributed species. In particular, introductions influence genetic diversity and population structure in the introduced range in unpredictable ways. Additionally, introgressive hybridization is often associated with successful establishment in new ranges. In this study, we explore the impact of introgressive hybridization on the polyploid Capsella bursa-pastoris in the New York City metropolitan area. We find Capsella bursa-pastoris in the New York City metropolitan area likely originated from multiple introductions from northern Eurasia, and that populations across the New York City metropolitan area are generally panmictic. As with Capsella bursa-pastoris in Eurasia, we discover evidence of introgression from the diploid Capsella rubella in this population. By evaluating ancestry in regions across the genome, we find introgressed regions are rich in gene content and contribute to genetic diversity in this population. These results suggest that introgressive hybridization before introductions may buffer species from the negative effects of population bottlenecks and allow for successful establishment.

Gene flow shapes the demographic stability and evolutionary potential of tropical forest trees, yet its dynamics may differ depending on the temporal scale at which it is assessed. We combined spatial genetic structure (SGS), parentage analyses, and reproductive success metrics to investigate historical and contemporary gene dispersal in four populations of Dicorynia guianensis across French Guiana, encompassing sites differing in environment and management history. A total of 1,528 individuals were genotyped using 66 nuclear and 23 plastid microsatellite markers, enabling high-resolution inference of biparental and maternal gene dispersal. Historical mating and dispersal parameters inferred from SGS revealed marked contrasts among populations. Some populations exhibited high historical gene dispersal distances and weak spatial genetic structure, whereas others showed stronger SGS and long-term aggregative dispersal patterns. Contemporary parentage analyses further highlighted differences in seed and pollen dispersal distances, parent assignment rates, and reproductive skew. In certain populations, pronounced reproductive inequality and reduced effective connectivity were observed, while others displayed more balanced reproductive contributions. By jointly evaluating long-term dispersal legacies and present-day reproductive patterns, our study demonstrates the value of combining indirect and direct genetic approaches to assess population dynamics and conservation status in tropical forest trees. This multi-temporal perspective provides a comprehensive basis for long-term monitoring and sustainable management in heterogeneous tropical landscapes.

Understanding relatedness in sharks is challenging due to uncertainty in distributions, low population densities and difficulties in sampling across life stages. In Fiji, bull sharks (Carcharhinus leucas), with an effective population size estimate of [~]258, aggregate at the Shark Reef Marine Reserve (SRMR), but gravid females disperse at the end of the year to give birth in adjacent rivers. Questions remain regarding reproductive connectivity, female returns across years, and kinship structure. Using population genomics on 296 bull sharks across age classes (neonates, young-of-the-year, juveniles, and adults) collected over a decade at the SRMR and in three adjacent rivers, we assessed familial connections. Direct genetic links, including first- and second-degree relationships, connected SRMR adults with young age classes in the Navua and Rewa rivers, providing evidence of reproductive connectivity. Within rivers, genetic similarities across cohorts revealed reproductive philopatry. Remarkably, several individuals sampled years apart were assigned to the same sire-dam pairs, indicating repeated pairings across breeding seasons. However, the few related links detected between the SRMR and the rivers may reflect incomplete sampling. Altogether, bull shark reproduction in Fiji seems influenced by reproductive philopatry and repeated pairings, suggesting added complexity in their reproductive behaviour.

Understanding how climate shapes intraspecific genetic turnover is critical for predicting biodiversity responses to global change, yet such analyses remain limited for systems where natural adaptation and human-mediated dispersal jointly structure diversity. Here, we investigate the spatio-temporal dynamics of genetic composition in the western honey bee (Apis mellifera) across Anatolia and Thrace, a major historical refugium harboring five subspecies. Using a dataset of 672 individuals genotyped at 30 microsatellite loci, we characterize population structure and model ancestry compositions as a function of environmental and geographic variables. We integrate Gradient Forests and Generalized Dissimilarity Modelling to identify key climatic drivers of intra-specific turnover and project future changes under multiple CMIP6 climate scenarios. We detect five major ancestral groups with widespread admixture structured by both spatial processes and environmental gradients. While geographic distance explains a substantial proportion of variation, climatic variables account for a large fraction of ancestry turnover. Spatial projections reveal distinct ecological regions corresponding to subspecies distributions, with high turnover zones aligned with major geographic and ecological barriers. Climate projections indicate substantial restructuring of ancestry compositions over the 21st century. Most ancestral groups show declines in persistence and resilience, whereas lineages associated with warmer and drier conditions expand under future scenarios. Regions of high uniqueness and refugia contract, while areas experiencing rapid turnover and novel ancestry compositions increase. Existing Genetic Conservation Areas provide incomplete representation of diversity and are projected to lose effectiveness under future climates. Our results demonstrate that climate change is likely to disrupt spatial genetic structure, promote admixture, and threaten persistence and resilience of honey bee populations. By modeling ancestry composition as a multidimensional proxy for genetic variation, for the first time to our knowledge, this study provides a scalable framework for forecasting intraspecific biodiversity dynamics and informing conservation and management strategies under global change.

The evolution of plant resistance naturally occurs in complex, multifaceted environments that consist of multiple simultaneous stressors. Understanding how shifting environmental contexts may shape resistance evolution requires empirical studies that consider the combined effects of interacting stressors on fitness and selection. Here, we examined how exposure to an insecticide impacts the evolution of resistance to the herbicide glyphosate in Ipomoea purpurea (common morning glory). Through a factorial field experiment, we manipulated glyphosate and an insecticide to estimate selection on glyphosate and herbivory resistance. We found that glyphosate acted as the primary agent of selection, favoring higher levels of glyphosate resistance. In the presence of glyphosate alone, positive correlational selection favored a combination of higher glyphosate and herbivory resistance, supporting prior work that suggested these traits may be linked. Importantly, insecticide exposure modified both glyphosate resistance and the strength of selection acting upon the trait by increasing resistance and weakening selection. Together, our results indicate that the evolution of herbicide resistance is context-dependent and that secondary stressors like insecticide can alter the evolutionary trajectories of plant defense.

Bacteriophages are not only the most ubiquitous biological entity on earth, they also display remarkable genetic diversity across and within populations. While macrodiversity has been extensively studied, the drivers of microdiversity (intraspecies genetic diversity) remain poorly understood, particularly in relation to phage lifestyle. The distinguishing ability of temperate phages to integrate themselves into the host genome has an unknown influence on the microdiversity present. This difference in microdiversity could impact the adaptability of phages to (a)biotic factors. To identify a possible association between microdiversity and lifestyle, we analysed 12 existing viromics datasets focusing on soil bacteriophages, including 41 412 viral genomes in total. We found that phages predicted to be temperate consistently exhibit significantly higher microdiversity than their virulent counterparts in eight of 12 datasets, whereas the remaining four datasets did not show a significant trend. The detected pattern holds across multiple quality thresholds and lifestyle prediction methods. These findings suggest that lysogeny may promote or preserve genetic variation within phage populations, with potential implications for phage-host coevolution and environmental adaptability.

Resource availability is a central driver of ecological and evolutionary processes, yet its effects on infectious disease and virulence are not fully understood. A key limitation is that many studies consider only a narrow range of resource conditions or a limited subset of host and pathogen traits, potentially obscuring non-linear relationships. Here, we quantify how a gradient of six food levels simultaneously shapes host fitness and pathogen performance in the Daphnia magna- Ordospora colligata system. Across two laboratory experiments, we measured infection rates, pathogen burden, host fecundity, survival, and filtration rates. Increased food availability enhanced pathogen fitness, with both infection rates and spore burden increasing with provisioning. In contrast, host responses were trait-specific: while fecundity increased with food availability, pathogen-induced reductions in fecundity (i.e., virulence) peaked at intermediate resource levels, despite continued increases in pathogen load. This pattern indicates that resource availability alters host tolerance as well as pathogen growth, generating non-linear disease outcomes. Host survival was unaffected by either food provisioning or infection, further demonstrating that resource availability can simultaneously influence host and pathogen traits in different directions. Our results highlight the importance of integrating multiple fitness components across provisioning levels to understand disease dynamics and suggest that ongoing anthropogenic changes in resource availability may alter host-pathogen interactions.

Understanding how viral communities vary across co-occurring hosts and environments is essential for assessing species-specific viral risks under changing land use and climate. This is particularly relevant for managing introduced bees, which face persistent viral threats themselves, as well as transmitting plant viruses. Here, we compare RNA viromes of the long-established honeybee (Apis mellifera, introduced to Tasmania in 1831) and the more recent invader, the bumblebee (Bombus terrestris, invasive since 1992), across 14 Tasmanian sites - an island still free of the viral vector, Varroa destructor. Using a metatranscriptomic approach on total RNA from whole bees, we identified insect- and plant-associated viruses and inferred phylogenetic patterns of insect viral sharing, divergence, and potential cross-species transmission. We also assessed spatial and environmental drivers of viral composition, diversity, and richness. Geographic longitude, precipitation, temperature, and pasture percentage influenced the total, insect-, and plant-associated viromes of B. terrestris. In contrast, for A. mellifera, only precipitation and temperature were associated with insect and plant viral alpha diversity and community composition. Phylogenetic analyses revealed that Black Queen Cell virus in A. mellifera from Tasmania has diverged from mainland Australian sequences, and two distinct sub-strains of Lake Sinai virus 1 were shared by both bee species. Lake Sinai virus 3 showed evidence of interspecies transmission between A. mellifera and B. terrestris. Notably, this study provides the first detection of Moku virus in Australian bees and globally in bumblebees, suggesting potential interspecies transmission among social Hymenoptera. Overall, our findings demonstrate local viral diversification and reveal that B. terrestris viromes are more strongly shaped by environmental factors than those of A. mellifera, underscoring the importance of monitoring invasive pollinators as reservoirs and vectors of viral emergence.

Extreme climatic events are reshaping ecosystems worldwide as individual organisms vary markedly in their ability to withstand these disturbances. Deciphering patterns of persistence on local scales is therefore critical for predicting biodiversity trajectories under intensifying climate extremes. In this study, we examined variation in thermal stress responses among individuals of the coral Stylophora pistillata species complex during a heatwave at Heron Island Reef, Australia. More than half of the focal coral colonies died on the reef, and survival of coral fragments maintained under ex situ common thermal stress conditions was significantly correlated with the survival of their source colony. This demonstrates that survival differences result largely from biological factors rather than differential thermal exposure across reef habitats. Under common garden conditions, we observed striking differences in bleaching severity and survival times among three sympatric cryptic taxa and their highly host-specific symbiont community. Within the most locally common taxon, corals from historically warmer and more seasonally variable reef habitats seem more susceptible to bleaching, contrary to expectations. Together, these results reveal how biological differences among cryptic taxa and among individuals can shape coral responses during a heatwave and advance our understanding of coral vulnerability in a rapidly warming world.

When invasive populations establish in regions far from their origin, they risk accumulating harmful mutations (genetic load) that limit population viability and subsequent spread. This may be exacerbated by the multiple, sequential bottlenecks experienced when invasions stem from a bridgehead population. However, populations may be able to purge genetic load when they can outbreed with other lineages from subsequent invasions. Here, we analyse global invasions of a species complex of persistently inbreeding ambrosia beetles, using genomic data (N=247) from invasive populations in Africa, North America and Australia, and from native populations in Asia. We focus particularly on one species of this complex (Euwallacea fornicatus) which poses a catastrophic threat to tree species worldwide and is rapidly expanding its global range. We uncover a single lineage of this species across South Africa, California and Western Australia, derived from an invasive bridgehead and containing almost no nuclear genetic variation. In South Africa we identify a second lineage that has repeatedly hybridised with the first lineage. Genetic patterns in the native range indicate that such opportunistic outbreeding may be common. Although purifying selection was evident in all lineages, native populations had fewer missense mutations than invasive populations, suggesting that opportunistic outbreeding may help purge fixed deleterious mutations when local lineage diversity is high. These findings show how inbreeding depression can affect populations even where inbreeding is common, and they highlight the biosecurity threat posed by subsequent gene flow into invasive populations.

Selective breeding in aquaculture is necessary to establish food security and meet demand for sustainably produced protein. An informed selective breeding program requires understanding how population structure, environmental adaptation, and human activities shape natural genetic variation in wild conspecifics. Unfortunately, wild variation remains poorly characterized for many commercially important aquaculture species. Here, we conduct the first range-wide study of genomic population structure for the eastern oyster (Crassostrea virginica) across thousands of miles (Texas, USA to Eastern Canada) using a 200K SNP array. We integrate population structure analyses, genotype-environmental associations, and structural variant detection to identify adaptive loci and quantify human-mediated genetic impacts. Our data confirms two ancestral clusters with a phylogeographic break between the Gulf and Atlantic (FST = 0.06) and highlights patterns of substructure within each region. We find evidence of unexpected patterns of genomic variation in two locations: evidence of Gulf ancestry in a mid-Atlantic estuary (Chesapeake Bay), and evidence of Atlantic ancestry in a Gulf estuary (Apalachicola Bay). While we cannot definitively determine the causes of these unexpected patterns, we show that they are consistent with direct and indirect human impacts in these estuaries. Genotype-environment association analyses with in situ temperature and salinity measurements were used to identify putatively adaptive loci, including SNPs within large structural variants (>1Mb). Our results identified genomic targets for aquaculture breeding programs aimed at climate resilience, reveal complex patterns of human impacts in managed systems, and demonstrate how seascape genomics can be used to improve aquaculture outcomes.